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The Large Eddy Simulation (LES) approach is adapted to the study of plasma microturbulence in a fully 
three-dimensional gyrokinetic system. Ion temperature gradient driven turbulence is studied with the GENE 
code for both a standard resolution and a reduced resolution with a model for the sub-grid scale turbulence. 
A simple dissipative model for representing the effect of the sub-grid scales on the resolved scales is proposed 
and tested. Once calibrated, the model appears to be able to reproduce most of the features of the free energy 
spectra for various values of the ion temperature gradient. 



I. INTRODUCTION 



Turbulence in plasmas shares several general features 
with fluid turbulence as modeled by the Navier- Stokes 
equation. In particular, microturbulence in a background 
magnetic field as described by the gyrokinetic formalism 
is thought to be characterized by a forward perpendic- 
ular cascade of free energy similar to the direct cascade 
of kinetic energy in the Richardson-Kolmogorov picture 
of fluid turbulence^El Although in plasma microturbu- 
lence, there does not exist an inertial range in the strict 
sense of the word, recent gyrokinetic simulations show 
that an asymptotically free, self-similar, and highly local 
cascade develops at high perpendicular wavenumbers 5 . 
While dissipative processes (due to a coupling to damped 
eigenmodes 6 ) are also active at low wavenumbers, a sig- 
nificant fraction of the free energy is transported to small 
spatial scales and dissipated there. In a simulation of 
the complete turbulent cascade process, it is thus im- 
portant, in principle, to capture all scales from the en- 
ergy injection range down to the smallest relevant dissi- 
pative scales. In certain situations, such direct numerical 
simulations (DNS) can be computationally expensive or 
even unfeasible, however. These DNS limitations have 
prompted the development of hybrid approaches mixing 
ab initio computation and modeling. In particular, Large 
Eddy Simulation (LES) techniques have been devised for 
simulating turbulent fluids at high Reynolds numbeiJ^. 
In these simulations, the large scales are computed explic- 
itly while the influence of the smallest scales is modeled. 
The aim of the present work is to extend this technique 
to the gyrokinetic equations which describe microturbu- 
lence in magnetically confined plasmas^^. 

As far as fluids are concerned, the main idea behind 
the development of LES techniques is the assumed exis- 
tence of a universal regime for the small scales. Indeed, 
the large scales in a turbulent fluid are very much influ- 
enced by the geometry of the flow. It is thus a priori not 
easy, and probably impossible, to design a general model 
for these large scales. On the contrary, if the Reynolds 
number is large enough, the smallest scales are supposed 
to be independent of the geometry and should only be 
affected by the physical properties of the fluid. Hence, 



there is a reasonable hope that the small scales can be 
represented by a general model. In practice, however, 
the Reynolds number is not always sufficiently large to 
reach such a regime and the geometry of wall-bounded 
flows has to be taken into account in most LES studies of 
turbulent fluids. Nevertheless, for almost half a century, 
LES have proven their ability to significantly decrease the 
numerical effort required to reproduce the main feature of 
large scale turbulent flows^l More recently, such meth- 
ods have also be en ap plied successfully to turbulence in 
conducting fluids^^. 

In the area of gyrokinetic turbulence, LES tech- 
niques have been explored, for instance, by Smith and 
HammettP^, considering hyperviscosity models for two- 
dimensional drift-wave turbulence - and the aim of the 
present study is to extend the LES methodology to gy- 
rokinetics in three spatial dimensions. In this context, it 
should be noted that plasma microturbulence is different 
from ordinary fluid turbulence in that it may be driven by 
various mechanisms such as, e.g., the presence of an ion 
temperature gradient (ITG) or an electron temperature 
gradient (ETG). It is thus to be expected that the mod- 
eling has to be adapted to the drive mechanism (rather 
than to the geometry of the system). In the present work, 
we will focus on the case of ITG turbulence. 

The paper is organized as follows. The fundamen- 
tal equations are discussed in Section [TTJ and the LES 



approach for gyrokinetics is presented in Section III to- 
gether with a simple dissipation model. The calibration 
of the model is discussed for Cyclone Base Case param- 
eters which is a standard ITG turbulence test case. An 
estimate of the truncated scales is proposed for various 



quantities in Section [IV] The robustness of the model 
in terms of parameter changes is analyzed in Section \V\ 



followed by a summarizing discussion in Section VI 



II. GYROKINETIC MODEL 

The LES approach can, of course, be studied in the 
context of a general gyrokinetic system, including mul- 
tiple particle species, electromagnetic fluctuations, colli- 
sions, general tokamak geometry, profile variations, and 
the like, as it is generally used in GENE! 14 " 16 [ How- 
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ever, in order to simplify the following discussions, we 
restrict to a reduced system here, working with only one 
ion species and treating the electrons as adiabatic. At the 
same time, a simple s-a model geometry is employed, and 
collisional effects are neglected. Moreover, the radially 
local version of GENE is used which s olves the gyroki- 
netic equations in a flux-tube geometr ^ 17 * 18 * , employing 
the field-aligned coordinates (x, y, z, vn , ji). The deriva- 
tion of the corresponding non-dimensional equations can 
be found in previous studies 16 . The resulting expressions 
read as follows. The time evolution of the ion distribution 
function fki in k space is given by 



dtfki = L[f ki ] + N[f ki , f ki ] + D[f ki ] 



(1) 



where the linear term can be split into three contribu- 
tions, L[f ki ] = L G [f ki ] + L c [fki] + L\\[fki], with 
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Here, hki is defined as the nonadiabatic part of the dis- 
tribution function, h ki = f ki + ZiFoifaToe/Toi where Zi 
denotes the charge number, F^ the background distribu- 
tion function, <j) the electrostatic potential and To e ,Toi 
the electron and ion temperatures. The first linear term 
Lq represents the influence of the fixed ion density (u n i) 
and temperature (uori) gradients, the second linear term 
Lc describes effects due to magnetic curvature, and the 
third linear term Ly contains the parallel dynamics in- 
volving magnetic trapping as well as the linear Landau 
damping. Meanwhile, the nonlinear term N represents 
the effect of the self-consistent electric field in the E x B 
drift of charged particles, 

N[fki, fki] = ^2 ( k x k v k * k y) J ^k'f{k-k')i , 
k' 

and the dissipation term D[f ki ] is given by 

D[f ki ] = ~ (a x (i k x ) n + a y (i k y ) n + a z d™ + a V]] <9™ ) f ki , 

where typically n = 4 is used, and the coefficients a x , 
a y , a z , and CLy,, can be adapted to each specific class of 
physical problems. In the local version of GENE used 
here, unknowns are Fourier transformed in the radial and 
poloidal directions, so that x and y are replaced, respec- 
tively, by k x and k y . The subscript 'fc' has been added to 
label such Fourier space quantities. Due to the imposed 
quasi-neutrality, the electrostatic potential (j) and the dis- 
tribution function are related via the linear equation 
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with A 2 = 2k\n/Bv and b { = v 2 , i k 2 _/(2Q 2 ci ). The func- 
tions Jo and To(bi) = exp(— bi) Io(bi) are, respectively, 
the Bessel and the scaled modified Bessel functions of 
order zero, where Q C i is the ion cyclotron pulsation and 
qi = Zie their charge. In the flux-tube geometry (sym- 
bolically defined by the metric coefficients 18 g xx , g xy and 
g vv \ the amplitude of the perpendicular wave vector k± 
is given by k\ = g xx k 2 + 2g xy k x k y + g yy k 2 and depends 
on z through the metric coefficients. Note that (<P)fs 
represents the flux surface average of the electric poten- 
tial. 

One property of the gyrokinetic equations which is of 
particular interest here is t he con servation of the free 
energy by the nonlinear terrrPE^Ql. rpj^ latter quantity 
is defined as 
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where h-ki — hi{—k Xl — k y ,z,v\\,t). The integration over 
the phase space of a given quadratic unknown \X k \ 2 = 
X(k x ,ky,z,v\\,fjb,t)X(-k x ,-k y ,z,v\\,fjb,t) is given by 



J dA k \x k \ 2 = ^^2Yl I ndzdv w d ^\ x > 



(4) 



where the sum over k™ s has to be understood as a 
sum from k x = (—N x /2 + 1) Ak x to k x = N x /2Ak x 
and the sum over ky NS corresponds to a sum from k y = 
(-N y /2 + l)Aky to ky = N y / '2 Ak y . Here, Ak x = 2tt/L x 
and Aky = 2i\jL y are the smallest wave vectors that can 
be used to represent periodic functions in rectangular do- 
main of size L x x L y . In practice, due to the symmetry of 
Fourier transform, negative k y modes are given by com- 
plex conjugation of positive k y modes. The volume V is 
defined in the chosen magnetic s — a equilibrium by 



(5) 
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The two-dimensional spectral density of the free energy 
is defined by 



h— k % f k i 
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(6) 



and the one-dimension spectral densities along k x or k y 
are simply given by 

—— ^ ^ £k x ,ky = ^ ^ £k x ,ky 

^DNS /c DNS 

The free energy balance can be expressed as follows 

d t £ = G-V. (7) 

The two terms in the right hand side represent the free 
energy injection Q and dissipation V. They are given by 



G = n 0i ^ 1 - [ dA k ^—^L G [f ki ], 

J-Oe J fOi 



(8) 



V = -n . 



^ / 'dk k ^D[f ki ). (9) 
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The free energy injection term Q is directly related to the 
ion heat flux, Qf. 



Q = u T i Qi 



(10) 



The heat diffusivity Xi an d nea ^ flux Qi = noiToiUriXi 
are considered as reference quantities for comparison be- 
tween gyrokinetic numerical solvers as well as with exper- 
iments. The appearance of heat flux as the free energy 
source stresses the importance of free energy balance in 
gyrokinetics. 



III. 



LES FOR GYROKINETICS 



The main objective of the LES technique is to ex- 
plore the large scale physics at a lower computational 
cost when compared to DNS. Reducing the cost of a gy- 
rokinetic simulation can be achieved by several ways. In 
a Eulerian approach, the distribution function is repre- 
sented on a fixed grid in five-dimensional phase space, 
using N x x N y x N z x N v ., x grid points. In a gyroki- 
netic LES, this grid is then to be coarsened. Considering 
that the perpendicular cascade processes are expected 
to transfer free energy from large spatial scales to small 
ones, we just employ the technique to the (x, y) grid in 
the present study. In the following, the DNS and LES 
grids correspond, respectively, to N x x N y = 128 x 64 
and N x x N y = 48 x 24, while N z 16, N V{] = 32, and 
= 8 are held constant. The coarsening of the grid can 
be viewed as a low pass filter, denoted hereafter by the 
— symbol, that sets to zero the highest k x and k y modes. 
The filtered distribution function will thus be labeled f ki . 
Applying this filter to the gyrokinetic equation yields: 

dtf ki = L[f ki ]+N[f ki J ki ]+D[f ki ]+T, (11) 

which contains a term T that depends explicitly on both 
the filtered distribution f ki and on the unfiltered distri- 
bution f ki 



T = N[f ki J ki ]-N[f ki J k i] 



(12) 



Except for the presence of T on the right hand side, 
Eq. (11) for f ki has the same form as Eq. |l]) for f ki . 
The term T is usually referred to as the sub-grid scale 
term, though in the present situation, the terminology 
sub-filter scale term would be more appropriate. In or- 
der to close Eq. ( 11 ), T must be approximated by a model 
that has to be expressed in terms of the filtered distribu- 
tion J ki : 



M[f k 



(13) 



The importance of the sub-grid scale term is illustrated 
in Fig. [I] Two free energy spectra are represented for the 
Cyclone Base Case 22 (CBC) for ITG driven turbulence 
(uj ni = 2.22, uj Ti = 6.92, q = 1.4, s = 0.796, e = 0.18, 
T e o/Tio = 1, Zi = 1). The perpendicular box sizes are 
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FIG. 1. Free energy: comparison between highly resolved 
DNS (black), and LES without model (blue) for the cyclone 
base test case. The resolution for the DNS is N x — 128 and 



N y = 64 and for the LES N x 



48 and N v = 24. 



given b y L x = L y = 125 pi. In GENE, numerical dis- 



sipation^ can be introduced via fourth-order derivatives 
along z and v\\: 



D[f k 



-dzdjfki 
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(14) 



where the values for the coefficients a z and a v ^ are to 
be adjusted appropriately. In Fig. [l] the DNS spec- 
trum is compared to the LES spectrum obtained by set- 
ting T = 0. Clearly, the free energy is piling up in the 
high k x range in the latter case due to the reduced high- 
wavenumber dissipation in the absence of small scales. As 
a secondary effect, the free energy appears to be pumped 
out more rapidly of the large scales where the gradient 
source term is active. Indeed, the free energy is trans- 
ferred to the small scales by the nonlinear term. These 
transfers have been identified as mostly local in Fourier 
spaced Consequently, the increase of activity of modes 
closer to the injection range can explain that the free 
energy is removed from the drive range more rapidly. 
Such a scenario is reminiscent of what is observed in un- 
derresolved DNS of Navier-Stokes turbulence, where the 
energy is also piling up in the large wavenumber range. 

In Fig. [I] the LES has been performed with T = 0, 
which can be considered as the simplest sub-grid scale 
model. Obviously, such a choice is too simple since the 
free energy spectrum deviates significantly from the DNS 
observations. The role of T can be understood by con- 
sidering the resolved free energy balance. Since f k i and 
f ki satisfy the same equation up to the term T, the free 
energy associated to f ki , referred to as the resolved free 
energy in the framework of a LES, must satisfy the fol- 
lowing equation: 



d t E = G-V-T 



(15) 



where the quantities £, 5, V are the same as £ , (5, D, 
except that they are defined using f ki and h k i instead of 
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FIG. 2. Sub grid contribution to the free energy balance com- 
pared with the free energy injection and dissipation terms. 
CBC parameters with a filter corresponding to N x = 48, 
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fki and hki- It should be noted, however, that all these 
global quantities are defined using a volume integration 
over dAk in which the sums are over fcjf s and ky ES and 
have to be understood as a sum from k x = (—N x /2 + 
1) Ak x to k x = N x /2Ak x and from k y = (-N y /2) Ak y 
to k y = Ny/2 Aky. Since the computational box sizes are 
the same in the LES and in the reference DNS, the same 
grid spacings Ak x and Ak y are used in both LES and 
DNS runs. However, the largest wave vectors are smaller 
in the LES than in the DNS: i^ ES = N x /2 Ak x < Kf s = 
N x /2Ak x and K^ s = N y /2Ak y < K™ 3 = N y /2Ak y . 
The new term 7^ is defined by 



Tt 



dAkUQi^^-^T (16) 

J-Oe ^Oi 



and represents the effect of the sub-grid scales on the 
resolved free energy. If the cascade picture applies, the 
effect of the should be to pump out the resolved free 
energy in order to mimic the transfer towards the unre- 
solved scales. If and hki are known from a DNS, the 
term T and consequently 7? can be computed exactly. 
Using the same parameter as in Fig. [lj has been com- 
puted and is shown in Fig. [2] It is indeed negative and 
represents a loss of resolved free energy. Its amplitude 
is compared to the resolved free energy injection rate Qj 
and dissipation rate Dj. On average, once turbulence is 
developed and a statistically stationary regime is reached, 
these three terms should be in balance, Qj ~ Vj + T^t. In 
the run corresponding to Fig.fTJ the ratio T^/Vj appears 
to be close to unity. Hence, Trie transfer of free energy 
between the resolved and the unresolved scales cannot be 
neglected. 

The development of models for representing the ef- 
fect of small, under-resolved scales on the large, resolved 
scales has been the subject of countless efforts in LES 
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FIG. 3. Resolved free energy spectra £ kx (top) and £ ky 
(bottom) obtained by varying the model coefficient compared 
with reference DNS. 



for fluid turbulence. However, the most commonly used 
models simply attempt to reproduce the transfer of ki- 
netic energy towards the unresolved scales by a dissipa- 
tive mechanism usually represented by an effective viscos- 
ity. Considering the analogy between fluid and plasma 
turbulence, it is proposed here to also use an effective 
dissipation which is modeled by the hyper- diffusion term 



M[f ki ] = -c ± k 4 ± h. 



ki • 



(17) 



It is easy to verify that such a model always gives a neg- 
ative contribution to the resolved free energy balance: 
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(18) 



The hyper- diffusion coefficient c±_ can be adjusted by 
comparing the results given by a reference (well resolved) 
DNS with results from an LES using Eq. (18). 
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In Fig. [3J the free energy spectra £ kx and £ ky , are dis- 
played for various values of c±. The black curve corre- 
sponds to the reference DNS run. It is compared to four 
LES runs corresponding to c± = {0,0.025,0.375,1.0}. 
Other values of c± have also been tested but are not 
shown here for clarity. All the spectra are time-averaged 
during the turbulent phase over a period of 2000Rq /vTi- 
Obviously, a too small parameter (c± = 0.025 dash- 
dotted) does not significantly improve the result when 
compared to the no- model case (c± = dotted). Also, 
a too large value of c± (x dotted) tends to over-damp 
the small scales which leads to an artificial accumulation 
of free energy in the large scales (small k). The optimal 
value appears to be close to c± = 0.375 (dashed). The 
corresponding LES reproduces fairly well the spectra of 
the resolved free energy both in k x and in k y . 



IV. ESTIMATE FOR THE SUB-GRID QUANTITIES 

The use of a model has been shown in the preceding 
section to improve significantly the agreement between 
DNS and LES in gyrokinetic simulations. However, since 
small scales are truncated in the LES runs, it is not pos- 
sible in LES to predict directly global quantities such as 
the total free energy, the total heat flux (or equivalently, 
the total free energy injection) and the total free energy 
dissipation. An estimate of the contribution from the 
truncated scales to these global quantities is certainly de- 
sirable if a comparison has to be made with experimental 
results. 

In this section, a simple estimate is proposed for the 
subgrid scale contribution to these quantities. It is noted 
that all these quantities (£, Q or P), generically repre- 
sented by Q, can be represented either by the their two- 
dimensional spectrum Q kx ^ k v or by their one-dimensional 
k x spectrum (Q kx ) and k y spectrum((3 fey ). In a LES, only 
the resolved part of Q, denoted hereafter Q is directly ac- 
cessible. It is given by: 

q= E Q kx = E Q hv - ( 19 ) 

\k x \<K^s \k y \<K^s 

The unresolved part of Q, denoted 5Q, contains three 
contributions 5Q = S X Q + S y Q + S xy Q: 

S *Q= E E Q kx,ky ' ( 2 °) 

\k x \>K%® \k y \<K%®i 

tyQ= E E Q kx,ky > ( 21 ) 

\k x \<K^ \k y \>K^ 

s xv q= E E Q kx,ky - ( 22 ) 

\k x \>K^\k y \>K^ 

In DNS, 5Q can be computed, but, in LES, it has to be 
estimated. Such an estimate can be obtained by noting 
that, in the large k x and k y ranges of LES runs, the 



quantity Q can often be approximated by decaying power 
laws: 



Qk x ^ A L,-ot x 
~ h x 5 



ky « A y k~ a y . (23) 



The amplitudes A x and A y as well as the exponents a x 
and oty can be estimated by linear regression from the 
LES spectra. In that case, the following estimates can 
be obtained : 
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(24) 
(25) 



Since these sums are finite, there is a priori no restric- 
tion on the values of the exponents. However, if the 
wave vector range is extended to infinity, these sums 
converge if and only if a x > 1 and a y > 1. Estimat- 
ing S xy Q is more difficult. However, assuming a sepa- 
rable spectrum Q kx ^ k v = q kx q 2 y ? it can be shown that 
5xyQ = S X Q SyQ/Q. In general, it is thus expected that 
the correction due to S xy Q is very small and can be ne- 
glected compared to 5 X Q or S y Q. This procedure has 
been used to estimate the total value of both the free- 
energy and the heat flux from the LES with the optimal 
value of c±_ and for the LES without model. 
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FIG. 4. Free energy injection spectra Q ky obtained with the 
model cj_ = 0.375, compared with reference DNS and without 
model. Unresolved spectra of LES c± = 0.375 and c± = 0.0 
have been c ons tructed using a regression method with the 
assumption (23). 



As far as the free energy is concerned, the best fit of 
the one-dimensional spectra yields a x = 2.21, a y = 1.74 
for c±_ = 0.375. These values are then used to com- 
pute the corrections S x £ and 5 y £. The LES estimate 
£les £ or total free energy can then be compared to 
the value measured from the DNS, £ DNS . It is found that 
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the LES estimate is in good agreement with the DNS 
value £ LES = 1.10 £ DNS . On the contrary, without a model 
(c_l = 0), the best fit of the free-energy one dimensional 
spectra yields a x — 0.79, a y = 0.85. In that case, the 
estimates p4|) and (25) would be divergent if the sums 



had to be extended to infinity. However, if the sums are 
limited to if DNS , it is possible to reconstruct the total free 
energy from the LES without model but the estimate is 
more than twice the value of the DNS: £ No MODEL = 2.1 £ DNS . 

The same procedure has been used for the free energy 
injection Q = ujti Qi- For c± — 0.375, the regression 
method yields a x = 3.60 and a y = 2.20, which gives 
the following estimate Q LES = 1.11 C/ DNS . Hence, again 
the value computed from the LES slightly overestimate 
the DNS value of the free energy injection. However this 
prediction is still in reasonable agreement with the DNS 
and provides a much better estimate than the no-model 
simulation for which £ NoMODEL = 1.38 £ DNS . 



V. ROBUSTNESS OF THE LES APPROACH 

The choice c± = 0.375 has proven to give a reasonable 
agreement between the LES and the DNS predictions, in 
the case of standard CBC parameters. However, the LES 
methodology is only useful if the model parameters don't 
have to be calibrated for each set of parameters. In this 
section, it is proposed to explore the robustness of the 
LES approach by varying the logarithmic temperature 
gradient ujtz- 

One of the most well known features of ITG turbu- 
lence is the Dimits shift 22 , i.e., a nonlinear upshift of the 
stability threshold with respect to a linear analysis. This 
upshift occurs when varying the values of the logarithmic 
temperature gradient ujti, while keeping the logarithmic 
density gradient uo n i constant. The explanation of such 
an effect is that turbulence nonlinearly transfers the free 
energy to the zonal flows (i.e., purely radial structures, 
corresponding to finite k x , but k y = fey = 0). These 
structures then suppress the ITG instability if its linear 
growth rate is not sufficiently large, and turbulence can 
not be driven even if the plasma is linearly unstable. 

In Fig. [5j all the parameters characterizing the CBC 
have been kept constant, except the logarithmic temper- 
ature gradient which is varied from the standard CBC 
value of uoti = 6.96 to ujt% — 5.5;6.0;8.0. All simulations 
are performed using the model described in Eq. (17), with 
c± = 0.375. The values uon = 6.0 and ujti = 5.5 are 
close to the nonlinear threshold. It is then observed that 
an important part of the free energy is stored into the 
zonal flows. Such a result is in qualitative agreement 
with the usual picture of the Dimits shifP^. On the con- 
trary when the temperature gradient is increased, the to- 
tal free energy increases, and peaks around k x pi ~ 0.05, 
kyPi ~ 0.15. 



Such a test shows that the LES approach can re- 
produce qualitatively the expected phenomenology at a 
much lower cost than the DNS. However, it is impor- 
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FIG. 5. Resolved free energy spectra £ kx (top) and £ ky 
(bottom) obtained with the model ( 17) cj_ = 0.375 for various 
values of the logarithmic temperature gradient uon- 



tant to assess the quantitative agreement between LES 
and DNS. For this reason, another comparison between 
LES and DNS has been performed for uj^i — 8.0. As 
shown in Fig. [6j the LES using the same model with 
c± = 0.375 again reproduces the resolved free energy 
spectra obtained from the DNS reasonably well. In par- 
ticular, there is a clear improvement when compared with 
the no-model run. The value of con — 8.0 corresponds 
to a more turbulent state than uj Ti = 6.96 and the LES 
appears to be fairly well robust in this regime. It should 
be acknowledged, however, that the situation is not fully 
satisfactory when turbulence strength is decreased. For 
instance, in Fig. [7] (ujtz = 6.0), the LES predictions, 
although still acceptable, starts to deviate significantly 
from the DNS results in the large scale range. This is 
reminiscent of a difficulty known in the development of 
model for LES in fluid turbulence. Very few models are 
capable to capture correctly the transition between lami- 
nar and turbulent flows. Probably, the very simple model 
proposed here also suffers from such a deficiency. 
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FIG. 6. Resolved free energy spectra £ kx (top) and £ ky 
(bottom) obtained with the model ( |17| ) c± — 0.375 for ujt% — 
8.0. Comparison with a highly resolved DNS with uti — 8.0 
and with LES without model. 



VI. DISCUSSION 

The study presented here shows that the concept of 
LES can be extended to three-dimensional (in space) gy- 
rokinetics. The very good agreement reported in Figs. [3] 
and [6] between fully resolved simulations and under- 
resolved simulations including a simple model for the fil- 
tered scales is encouraging. It shows that the model, 
calibrated for a given value of the temperature gradient 
in ITG turbulence, is able to reproduce the large scale 
spectra of the free energy for higher temperature gradi- 
ent. It should be noted that such a test is quite demand- 
ing since the free energy has to be reproduced for each 
scale. Hence, the model has not only to be able to dis- 
sipate the correct amount of free energy, it also has to 
distribute the dissipation correctly amongst the different 
scales. 

It should be acknowledged, however, that the robust- 
ness of the model has not been proved in the most general 
sense. For instance, when the turbulence level is lowered 



FIG. 7. Same as Figure [6] with ujt% = 6.0. 

by decreasing the parameter uoti, the agreement between 
the under-resolved and the fully resolved simulations be- 
come less and less satisfactory - although such behavior 
can be understood and even anticipated since, for low 
co>Ti, turbulence is not fully developed and the cascade 
picture starts to break down. An interesting extension of 
the present approach would then be to apply the dynamic 
procedure used to calibrate automatically the amplitude 
of sub-grid scale models in LES for fluid turbulence^. 
The dynamic procedure is known to be able to predict 
the transition between turbulence and laminar flows by 
automatically setting the model amplitude to zero in the 
laminar regime in which the small scales are not active. 
Hence, although the models are designed by using con- 
cepts valid for fully developed turbulence, the dynamic 
procedure seems to be able to extend their validity into 
totally different regimes. However, the implementation of 
the dynamic procedure in gyrokinetics is more intricate 
than the use of the simple model studied here. 

Finally, it is interesting to discuss the computational 
gain obtained in the LES simulation presented in the pre- 
vious section. Although for the present case, the fully re- 
solved reference simulation is not using a very large grid, 



8 



the LES grid (and hence the required memory) can be 
reduced by 86%. In terms of CPU time, the gain is even 
higher. Indeed, the simulation can be performed with a 
larger time step since the smallest scales are larger than 
in the reference simulation (in practice, with the grid res- 
olution chosen in this study, the time step is increased by 
a factor of about two in the LES when compared to the 
DNS). As a consequence, the overall computational effort 
required for the LES runs appears to be more than 20 
times smaller than in the DNS simulations. This finding 
indicates that the LES approach is quite promising in the 
context of gyrokinetics. Further studies along these lines, 
also exploring other types of plasma turbulence (driven, 
e.g., by trapped electron modes or electron temperature 
gradient modes), will be the focus of future work. 
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